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. Fresnel and Kirchhoff made the proofs for Huygens' principle for the prop- 

D , agation of light respectively. In this paper, another proof for this principle 

^ I is given, in which, by using Poisson's formula for the three-dimensional wave 

CN I equation with initial-value conditions, the back-travelling waves are cancelled 
by the conterterms. 

■ PACS. 4210-Optical propagation and transmission in homogeneous media 
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On I 1. Introduction 

^ ' According to Huygens' principle for the propagation of light, every point 

O ■ illuminated by an incident primary wavefront becomes the source of a sec- 

^ . ondary wavelet such that the primary wavefront at a later time results from 

the superposition of amplitudes of all wavelets. It is generally agreed that 
1^ I back-travelling waves would appear from the above given principle. To dis- 

' prove this, Fresnel and Kirchhoff made a proof for the principle[l, 2]. Now 

we try to give another proof for the principle that when a point light source 
at Po, emits a series of spherical wavelets, the effect at a point P could be 
^ i found in either of two ways. First allow one of these wavefronts to proceed 

^ [ until it reaches P. Second, divide one wavefront into small areas, so that each 

area can generate a wavelet. Each of these wavelets then produces effects at 
P which can be summed, according to the principle of superposition, to give 
the total effect of the wavelets at P. These two methods must yield identical 
results. For this purpose, we should state Huygens' principle as follows. 
Suppose that the initial-value problem 
Lu{x, y,z,t) = t > 0, — oo < x,y, z < oo 
u{x,y,z,0) = (p{x,y,z,) {x,y,z)eD 
ut{x,y,z,0) = ^p{x,y,z,) {x,y,z)eD 
where L is an operator and D is a domain in R^, has the unique solution 
u{x, y, z, t) = f{x, y, z, t) {x, y, z) E Dt, 
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Here is the domain in at time t which is enclosed between the front 
and rear edges of the wave. Taking ti > as the initial instant of time, then 
the unique solution of the following initial-value problem 

Lu{x, y, z,t) = 0, t > ti, —oo < x,y, z < oo 

u{x, y, z, ti) = f{x, y, z, h) (x, y, z) G Ai 

ut{x, y, z, ti) = ft{x, y, z, ti) {x, y, z) e Ai 
is also f{x, y, z, t) {x, y, z) e Df. In particular, if we take t — t2{> ^i), then 

u{x, y, z, t2) = f{x, y, z, ta) {x, y, z) G A2 

That is to say, any point on a wave surface at ti behaves like a new source 
from which secondary waves spread out in a way satisfying the original wave 
equation. 

To make this obvious, we first analyse the one-dimensional wave equation 
with initial- value conditions. 

2. Initial- value problem of the one-dimensional wave equation 
Consider the initial-value problem of the one- dimensional wave equation 



Utt = a^Uxx t > 0, —00 < X < 00 

u{x, 0) = ip{x) — 00 < X < cxo (1) 

Ut{x, 0) = ip{x) — 00 < X < 00. 

The solution of the problem, by D'Alembert's formula, is 

1 1 r^+o-i 

u{x, t) = -[ip{x + at) + ip{x - at)] + — J ^ i/j{Od^ (2) 

To simplify the problem, let ip{x)=0. Then, at t=ti, we have 

u{x, ti) = ^[(fix + ati) + (p{x - ati)]. (3) 

The function (p{x — ati) is a wave moving to the right with a speed a and 
the function (p{x + ati) is a wave moving to the left with a speed a. 

The initial perturbation (p{x.) for ti > splits, so to speak, into two similar 
perturbations {^)<p{x ± ati), each with half of the intensity. (see Fig. 1,2) 
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Fig. 1 



Fig. 2 



Consider the case when t = t2 > ti again. On the one hand, by formula 
(121), we have 



On the other hand, we may take u{x, ti) as initial perturbation, and form 
the initial-value problem 



where f{x, ti) = ^[ip{x + ati) + {p{x — ati)]. 

We transform the initial- value problem (5) to an initial-value problem in 
which the initial time is r = 0. To this end we will introduce a new variable 
r = t — ti instead of the variable t. Under this change of variable, the 
initial-value problem (5) for the function 



u{x, = - [ip{x + at2) + ip{x - at2)] 



(4) 



Utt = a^Uxx t > ti, — oo < X < oo 
u{x,ti) = f{x,ti) — oo < X < oo 

Ut{x,ti) = ft{x,ti) — oo < X < CO, 



(5) 



u{x, t) = u{x, T + ti) 



becomes 



'TT 



•XX 



r > 0, — oo < x < oo 



u{x,0) = f{x,ti) 

Urix,0) = ftix,ti) 



— oo < X < oo 



— oo < X < oo 



(6) 



Using D'Alembert's formula again, we have 
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Let t = t2. Then 



U{x,t2) = U{x,t2 - ti) 

1 1 

-f[x + a{t2 - + -f[x - a{t2 - 



+ 



2' 
1 

2a 



x+a{t2—t-i) 



' x—a{t2—ti) 

Substituting f{x, ti) = h[ip{x + ati) + ip{x — ati)], we find that 



(7) 



u{x, t2) = -^fi^ — cih) + '^fi^ ~ 2ati + at2) + -^^{x + 2ati — 0^2) 
11 1 

+ -V^lx - at2) - ^V^l^; + 2ati - at2) 



(8) 



We can easily see that the initial perturbations f{x,ti) = {l/2)[ip{x + 
atl) + y^ix — atl)] for t2 > ti split into eight perturbations. The first four 
terms rise from the initial value f{x,ti), (as shown in Fig.S.a). The last four 
terms rise from the initial speed ft{x, ti).{Fig.3.h) 



j!p{x + at2) + 2ati — at2) ^'fi(x — 2ati + at2) ^ip(x — at2) 







~at2 —atl a{t2 — 2ti) a(2ti — t2) ati at2 

Fig. 3. a 

Formula (jH)), and Fig. 3. a , and Fig.S.b show there are back-traveling waves 
in the secondary wavelets, i.e. \ip{x + 2ati —0^2) and |v9(a; — 2ati + 0^2), rising 
from the initial value f(x,ti). However, there are also two counterterms, 
— — 2ati + 0^2) and —\ip{x + 2ati — 0^2), rising from the initial speed 
f((x,ti). So the back-traveling waves are cancelled by each other. The final 
result is 



u{x, 12) = ^ [ip{x + at2) + (f{x - at2)\ 



(9) 
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^ip(x + at2) 



\ip(x - 0*2) 




Fig. 3.C 



which agrees with formula (4). (as shown in Fig.3.c) 
The proof is similar when ^(x) 7^ 0. 

3. Proof of Huygens' principle 

Let Po(Fig.4) be a point light source at which spherical monochromatic 
scalar waves are being generated continually beginning at t=0. At t = ti > 0, 
the sphere of radius cti, with centre at Pq, and called S^°, is full of light, 
where c is the velocity of hght. The light disturbance at a point Q within 
the sphere or on the spherical surface may be represented by 
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„, ^ , Asinicuti - kr) 

f[r,ti) = , 

r 

where A is the amphtude at unit distance from the source, and r = |PoQ|- 
As usual, cu is the angular frequency and k — cu/c. 

We shall determine the light disturbance at a point P at ^2 > ^i- On the 
one hand, by allowing the light to proceed until it reaches P, the disturbance 
must be 

= (10) 

where R = \ PqP \ . 

On the other hand, the "light sphere" S*^" may be regarded as the initial 
disturbance and the scalar wave satisfies the wave equation □/ = 0, where 

□ - — -c^f— — —) 
is the D'Alembert's operator. Therefore we form the initial-value problem 

□it(x, y, = t > ti, —oo<x,y,z<oo 
u{x,y,zM)^ f{rM) {x,y,z)eS^° (11) 
ut(x, y, z, h) = ft(r, h) (x, y, z) e S^^ 

Let T = t — ti and u{x, y, z, t) = u{x, y,z,T + ti) = u{x, y, z, r). Then the 
above initial- value problem becomes 

Du{x, y, z.t) = r > 0, —oo < x,y, z < oo 
u{x,y,z,0)^f{r,h) {x,y,z)eSS', (12) 
Ur(x, y, z, 0) = ft{r, h) {x, y, z) e 

According to Poisson's formula for the three-dimensional wave equation 
with initial- value conditions [3, 4], the above initial- value problem has the so- 
lution 

^iP^r) = ^[^ J I kMds+ J I h{r,h)ds] (13) 

qP cP 

where S^^ is a spherical surface, centered on P, with a radius p=cr.(Fig. 4) 
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To evaluate the two integrals on the right-hand side of formula (jl3|) . we 
construct a sphere of radius r(r < cti) with center Pq . This sphere cuts the 
sphere S^^ in a circle c centered on PqP, and passing through a point Q on 
the spherical surface 5*^. Let the light disturbance at the arbitrary point on 
the circle c be represented by f{r,ti) = Asin{ujti — kr)/r. Another circle c' 
can be constructed using a sphere, centered on Pq, of radius r+dr. This circle 
is also centered on PqP and passes through a point Q' on the S*^. There is a 
ring zone between the c and the c' on the S"^. Let the area of the ring zone 
be ds, and let (p be the angle between PqP and PQ. From Fig. 4 

r2 = p2 ^ i?2 _ 2pRcosip. 

So that 

rdr = pRsimpdip 

and therefore 

ds = 27ip sin (pdip 

Now we must discuss two cases. 
Case I R + p < cti. 



zirpr 
R 



dr. 



(14) 
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Then the two integrals in right-hand side of formula (jl3p are taken over 
the whole spherical surface S^^ and r varies from R — ct to R + ct. It follows 
that 



u[p,Tj = — — / - ■ — simujti — krjZTT—ar 

A-KC OT Jn^cT P R 

+ / - ■ cos(co'ti — krjZTT — ar\ 

Jr-ct P r R 

pR+CT pR+CT 

I sin(a;ti — kr)dr + / u;cos(a;ti — kr)dr\ 

JR-CT JR-CT 



A , d 



2cR'dT 

= — [sin(ct;ti — kR — kcr) + sin(co'ti — kR + kcr) 
2R 

— sin(c(Ji(:i — kR — kcr) + sin(co'ti — kR + kcr)]. 
Substituting kc = u and let t = t2, we find that 



''^(P) ^2) = — [sin(2co'ti — ujt2 — kR) + sin(co't2 "~ kR) 
2R 

— sin(2u;i(:i — ujt2 — kR) + sin(u;t2 — kR)] (15) 

= — sin(co't2 — kR). 
R 

Case 11 R + p > cti. 

Then the two integrals in the right-hand side of formula (jl3p are taken 
over a portion of the spherical surface S^and r varies from R—ct to R+ct—'j, 
where < 7 < 2cr and i? + cr — 7 = cti. It follows that 

A d r^+(^^-'y 
u{p,t) = ^^[^y sm{ujti - kr)dr 

pR+CT — f 

+ / iu cos{ujti — kr)dr] 

JR-CT 

= — [sin(co'ti — kR — kcr + k'y) + sin(co'i(:i — kR + kcr) 
2R 

— sin(a;ti — kR — kcr + k'j) + sin(ci;ti — kR + kcr)]. 
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Substituting kc = uo and let t = t2, we find that 



u{p, ^2) = 777^[sin(2co'ti — U!t2 — kR + k^y) + sin(ct;t2 — kR) 
2R 

— sin(2ci;ti — u;t2 — kR + k'j) + sin(ti;t2 — kR)] (16) 

= — sin(ct;t2 — kR). 
R 

Both formulae (fT3j) and ()16|) agree with the formula (fTUj) . 
We can rewrite the two terms in the right-hand side of the first equal-sign 
in formula (jl5|) . Thus 

— sin(2u;ti — ujt2 — kR) 
2R 



and 



become 



and 



-—— sin(2c<jti — ut2 — kR) 
2R 



-— sinfc[i? + c(t2-2ti)] 
Zri 

— smk[R + c{t2-2ti)] 



respectively. We can see that the terms are the back-traveling waves and 
they cancel each other. The case is similar in the formula ()16|). 

It should be note that, more generally if we let the primary wave be 

/(r - ct) 
r 

where / is an arbitrary function, we can also obtain the same result. 
13-mail: zxliu2002@hotmail.com 
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